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ABSTRACT 



We present Kepler observations of the bright (V=8.3), oscillating star HD 
179070. The observations show transit-like events which reveal that the star is 
orbited every 2.8 days by a small, 1.6 -REarth object. Seismic studies of HD 179070 
using short cadence Kepler observations show that HD 179070 has a frequency- 
power spectrum consistent with solar-like oscillations that are acoustic p-modes. 
Asteroseismic analysis provides robust values for the mass and radius of HD 
179070, 1.34±0.06 M© and 1.86±0.04 R© respectively, as well as yielding an age 
of 2.84±0.34 Gyr for this F5 subgiant. Together with ground-based follow-up ob- 
servations, analysis of the Kepler light curves and image data, and blend scenario 
models, we conservatively show at the >99.7% confidence level (3(7) that the tran- 
sit event is caused by a 1.64±0.04 REarth exoplanet in a 2. 785755±0. 000032 day 
orbit. The exoplanet is only 0.04 AU away from the star and our spectroscopic 
observations provide an upper limit to its mass of ~10 MEarth (2-c). HD 179070 
is the brightest exoplanet host star yet discovered by Kepler. 

Subject headings: planetary systems — stars: individual (HD 179070: KIC3632418) 
— stars: oscillations — stars: interiors — stars: late-type - stars: activity - stars: 
magnetic field — techniques: photometric — Facihties: The Kepler Mission 
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1. Introduction 

The NASA i^ep/er mission, described in Borucki et al. (2010a), surveys a large area of 
the sky spanning the boundary of the constellations Cygnus and Lyra. The prime mission 
goal is to detect transits by exoplanets as their orbits allow them to pass in front of their 
parent star as viewed from the Earth. Exoplanets, both large and small and orbiting bright 
and faint stars, have already been detected by the Kepler mission (see Borucki et al. 2010b; 
Koch et al. 2010a; Batalha et al. 2010; Latham et al. 2010; Dunham et al. 2010; Jenkins et 
al. 2010a; Holman et al. 2010; Torres et al. 2010; Howell et al. 2010; Batalha et al. 2010; 
Lissauer et al. 2011) and many more exoplanet candidates are known (Borucki et al. 2011). 

We report herein on one of the brightest targets in the Kepler field, HD 179070 (KIC 
3632418, KOI 975). This V=8.3 F6IV star has been observed prior to the Kepler mission 
as a matter of course for essentially all of the Henry Draper catalogue stars. Catalogue 
informatioiij^ for HD179070 provides a T^s of 6137K (spectral type F6 IV), an Hipparcos 
distance of 108±10 pc, [Fe/H]=-0.15, an age of 2.8 GYr, a radial velocity of -28 km/sec and 
E(b-y)=0.011 mag. The age, metallicity, and color excess reported in the catalogue come 
from the photometric study by Nordstrom, et al. (2004). nothing special to distinguish it 
from any other stars in the HD catalogue. 

HD 179070 is strongly saturated in the i^ep/er observations. Techniques to obtain good 
photometry from saturated stars are employed by the Kepler project and have been used to 
analyze light curves of other Kepler saturated stars (e.g., Welsh et al. 2011) While a star 
will saturate the CCD detectors on Kepler if brighter than V~R~11.5, collecting all the 
pixels which contain the starlight, including the bleed trail, allow a photometric study to 
be performed (see Gilliland et al. 2010a). Kepler light curves of HD179070 from the early 
data as well as subsequent quarters of observation have shown a repeatable transit-like event 
in the light curve. We will examine the nature of this periodic signal using Quarter to 
Quarter 5 Kepler observations and will show that it is caused by a small exoplanet transiting 
HD 179070 with a period near 2.8 days. We will refer to the host star as HD179070 and the 
exoplanet as Kepler-21b throughout this paper. 

Our analysis for this exoplanet follows the very complete, tortuous validation path fully 
described in Batalha et al. (2010). The interested reader is referred to that paper for the 



University of California and the California Institute of Technology, the Mayall telescope at Kitt Peak National 
Observatory, and the WIYN Observatory which is a joint facility of NOAO, University of Wisconsin-Madison, 
Indiana University, and Yale University. 

^http:/ /simbad.u-strasbg.fr/simbad/sim-fid 
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detailed procedures which we do not completely repeat herein. Section 2 will discuss the 
ifejiZer observations including both the light curves and the pixel image data. §3 & §4 discuss 
ground-based speckle and spectroscopic follow-up observations we have obtained and Section 
5 presents the asteroseismic results. Using a full analysis of the observations, we produce a 
transit model fit for the exoplanet, determine many observational properties for it, and end 
with a discussion of our results. 



2. Kepler Observations 

2.1. Kepler Photometry 

The Kepler mission and its photometric performance since launch are described in 
Borucki et al. (2010a) while the CCD imager on-board Kepler is described in Koch et al. 
(2010b) and van Cleve (2008). The Kepler observations of HD 179070 used herein consist 
of data covering a time period of 460 days or Kepler observation Quarters through 5 (JD 
2454955 to 2455365). The photometric observations were reduced using the Kepler mission 
data pipeline (Jenkins et al. 2010b) and then passed through various consistency checks 
and exoplanet transit detection software as described in Van Cleve (2009) and Batalha et 
al. (2011). Details of the Kepler light curves and the transit model fitting procedures can 
be found in Howell et al. (2010), Batalha et al. (2011) and Rowe et al. (2006). 

The top panel of Fig. 1 shows the raw Kepler light curve of HD 179070 where the larger, 
low frequency modulations do not likely represent real changes in the star. Thermal jumps 
in the focal plane temperature near day 100, 120, and 370 are apparent (see Van Cleve 2009) 
due to safeing events in Q2 and Q5. Normalized and phase-folded Ught curves are produced 
and the transit event in the phased light curve is modeled in an effort to understand the 
transiting object. The middle panel of Fig. 1 shows a blow up (as highlighted in yellow) of 
a section of the raw light curve and represents a typical normalized result. The dotted lines 
mark the individual transit events and the entire Kepler light curve shown here contains 
164 individual transits. Taking random quarters of the total light curve and binning them 
on the transit period provide consistent results as shown in the bottom panel of Fig. 1 
albeit of lower S/N. The bottom panel (Fig. 1) shows the entire phase folded light curve 
after detrending and binning all available data. Each bin has a width of 30 minutes and 
the red curve shows our transit model fit to the data. Points marked with 'o's show where 
the exoplanet occultation would occur for a circular orbit (i.e., light curve phased at 0.5) 
or where evidence of a secondary eclipse would be seen if the event arises from a blended, 
false-positive eclipsing binary. 
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Due to the sparse (30-minute) sampling of the image data over the 3.4 hour long transit 
as well as the low S/N of any individual transit event, a search for transit timing varia- 
tions produced no detectable periodic signal with an amplitude of 8 mins or greater and no 
significant deviations at all for any period. 

2.2. Kepler Images and False Positive Analysis 

HD 179070 is saturated in Kepler images. In a typical quarter the photometric aperture 
covers 128 pixels, of which 39 are typically saturated in each image. A direct pixel image of 
HD 179070 in quarter 5 is shown in Figure |2j 

False positive identification for unsaturated Kepler targets proceeds by forming an aver- 
age difference image per quarter by subtracting an average of the in-transit pixel values from 
average nearby out-of-transit pixel values. The resulting difference image is used to compute 
a high-accuracy centroid of that difference image (Torres et al. 2010). For unsaturated tar- 
gets the difference image provides a star-like image at the actual location of the transiting 
object. The high- accuracy centroid is computed by performing a fit of the difference image 
to the Kepler pixel response function (PRF; see Bryson et al. 2010). When the difference 
images are well behaved, the centroid can have precisions on the order the PSF scale divided 
by the photometric signal to noise ratio of the transit signal. For Kepler-21b the SNR per 
quarter of about 25 would support centroiding precisions for sources near the primary of 
about 0.06 pixels. Having lost spatial resolution due to saturation the centroiding capability 
will be degraded compared to this. 

The high degree of saturation of HD 179070 prevents the direct application of the above 
centroiding technique, but we believe that a PRF fitted centroid to the non-saturated pixels 
in HD 179070's wings provide useful, albeit less accurate, results. 

Figure [3] shows the average difference image for Quarter 5, which is typical of other 
quarters. The pixel values which change during a transit are at the ends of the saturated 
columns as the amount of saturation that spills out of these ends decreases during transit. 
We also detect the (in phase) transit signal in the stellar wings around the core to the right 
of the saturated columns. Some saturated pixels become brighter on average during the 
transits which we ascribe to negative difference values in the pixel-level systematics linked 
to the non-linear behavior of saturated pixels and to large outliers caused by image motion 
events. Such negative difference image pixel values are commonly associated with saturated 
Kepler targets in which we observe shallow transits. 

We can enhance our view of the stellar wings around the saturated core of HD 179070 
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by setting pixels in columns that have saturated pixels in the direct image to zero in the 
left panel of Figure |4} The right panel shows the modeled difference pixel image created by 
simulating the transit on HD 179070 using the PRF model of Bryson et al. (2010) and the 
measured transit depth, similarly setting the saturated columns to zero. We see that for the 
stellar wings on the right side there is a reasonable qualitative match between the modeled 
difference image and the observed pixel values. However, on the left side, the observed wings 
have slightly negative values. Since we can not know the exact center of the stellar image 
due to the many saturated pixels, we expect that the PRF-fitted centroid will have some 
x,y bias. These negative values appear for all quarters, but their locations around the core 
vary significantly from quarter to quarter, indicating that the PRF-fitting bias varies from 
quarter to quarter. 

The PRF fit is performed via Levenberg-Marquardt minimization of the difference 
between the modeled and observed difference values, using non-saturated pixels in the model 
image whose value exceeds 10~^ of the summed model pixel values. This fit is done on both 
the difference and direct image, allowing us to compare the centroid of the difference image 
with the measured centroid of HD 179070 using the same pixels. The resulting in-transit 
pixel offsets are given in Table |2| with their formal uncertainties, for quarters 1, 3, 4 and 
5 (no quarter 2 data are available). These uncertainties do not include the expected but 
unknown fit bias due to the negative values in the observed difference image. In quarter 1 
we see centroid offsets exceeding a pixel, but this is almost certainly due to significant image- 
motion related to pixel-level systematics that were eliminated in later quarters (Jenkins et 
al. 2010c). In the remaining quarters the offsets are smaller, particularly in quarters 4 and 
5 where the offsets are less than a pixel, and there is no consistency, as expected, in the 
offsets. 

To gain confidence in PRF fitting using pixels in the stellar wings and in using the Kepler 
image data for HD 179070 as a way to set limits on possible background sources which could 
make this a false positive, a series of modeled centroids was produced. The source of the false 
positive event was modeled to be caused by a very dim variable test star (e.g., a background 
eclipsing binary) with a Kepler magnitude of 18 in which a variation (eclipse) was added 
with a depth of 50%. Models were produced with the faint test star placed at various 
pixel positions near the expected center of HD 179070 but such that some of the test stars' 
light would spill into the unsaturated wings. While the models did not exactly reproduce 
the observed data due to unknown systematics, the location of the saturated columns was 
faithfully reproduced when the transit was assumed to be on HD 179070. When the transit 
was placed on a background star offset by more than half a pixel in column the saturated 
column clearly became inconsistent with the data. When the transiting object is offset 
by a pixel or more in the row direction the location of the wings in the model data became 
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similarly inconsistent with the observed data. In this way the models show that the observed 
difference image can clearly rule out the possibility that the transiting object is located more 
than a pixel (4 arcsec) from HD 179070. Row offsets displace the model stellar wings in a 
difference image by large fractions of a pixel relative to the observed difference images. The 
effect of a column offset is shown in Figure [5} where we see that a column offset of one pixel 
(or more) causes the transit signal in the difference image to disappear from one or more 
of the saturated columns again not consistent with the Kepler observations. Based on our 
PRF model results, we are confident that the source of the transiting object is <1 pixel of 
the x,y position of HD 179070. 



3. High Resolution Imaging 

3.1. Speckle Observations 

A major part of the Kepler follow-up program (Batalha et al. 2010) used to find false 
positives as well as provide "third light" information to aid in Kepler image analysis is speckle 
imaging. We perform our speckle observations at the 3.5-m WIYN telescope located on Kitt 
Peak where we make use of the Differential Speckle Survey Instrument, a recently upgraded 
speckle camera described in Horch et al. (2010). Our speckle camera provides simultaneous 
observations in two filters by employing a dichroic beam splitter and two identical EMCCDs 
as the imagers. We generally observe simultaneously in "V" and "R" bandpasses where "V" 
has a central wavelength of 5620A , and "R" has a central wavelength of 6920A , and each 
filter has a FWHM=400A. The details of how we obtain, reduce, and analyze the speckle 
results and specifics about how they are used eliminate false positives and aid in transit 
detection are described in Torres et al. (2010), Horch et al. (2010), and Howell et al. (2011). 

The speckle observations of HD 179070 were obtained on 17 Sept. 2010 UT and consisted 
of one set of 1000, 40 msec speckle images. Our R-band reconstructed image is shown in 
Figure 6 and along with a nearly identical V-band reconstructed image reveals no companion 
star near HD 179070 within the annulus from 0.05 to 1.8 arcsec to a limit of {5a) 5.3 
magnitudes fainter than the target star (that is brighter than R~13.6). At the distance of 
HD 179070 (d~108 pc), this annulus corresponds to distances of 5.4 to 194 AU from the 
star. We note that any stellar companions or massive exoplanets (~20 Earth masses or more) 
inside 5.4 AU would be easily detectable in the radial velocity signature (see §4). We note 
that while reaching to 5 magnitudes fainter than the target star eliminates bright companions 
and some fraction of low-mass faint associated companions, it does not completely rule out 
the probable larger population of faint background EEs (see §7). 
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3.2. MMT AO Imaging 

AO images of HD 179070 were obtained on 2010 September 23 UT using the ARIES 
instrument on the 6.5-m MMT. ARIES is a near- infrared diffraction-limited imager, and was 
operated in the F/30 mode (0.02" per pixel) in both the J and Ks filters. The combined Ks 
image was created by combining 18 images: one 1-s and 17 0.9-s exposures, with two initial 
pointings at each exposure time and a raster of 4 4-point dithers with a jittered 2" offset 
between each position. The combined J image was similarly made from 17 0.9-s exposures 
(one single and sixteen dithered images). The images were shifted, sky-subtracted, and 
combined using xmosaic in the IRAF package xdimsum. The final AO images are shown in 
Figure 7. 

The seeing was relatively poor, with image FWHM of 0.33" in J and 0.29" in Ks. A 
faint companion is detectable in the Ks data, and hinted at in the J data. PSF fitting, using 
an analytical Gaussian model, found a magnitude difference of A — Ks = 3.7 + / — 0.1, while 
a lower limit was estimated for A — J > 4.0. The separation between the two components 
was found to be 0.7" +/- 0.05", with an approximate position angle of 135 degrees east of 
north. 

The faint companion was detected right at the limit of the ARIES frame, and at a 
distant of just over 2 FWHM. At the object's distance, 0.7", the estimated Ks detection 
limit for additional companions was 4.2 mag fainter (3.6 mag in J), increasing to 5.7 mag in 
Ks (5.1 in J) at 1", 7.5 mag in Ks (7.2 in J) at 2", and 8 mag in Ks (7.8 in J) beyond 4". 
In order to get better magnitude limits and a firm J-band detection, additional AO images 
were acquired using Keck. 

3.3. Keck AO Imaging 

Near-infrared adaptive optics imaging of HD 179070 was obtained on the night of 22 
February 2011 and 23 February 2011 with the Keck-II telescope and the NIRC2 near-infrared 
camera behind the natural guide star adaptive optics system. NIRC2, a 1024 x 1024 HgCdTe 
infrared array, was utilized in 9.9 mas/pixel mode yielding a field of view of ~ 10". Observa- 
tions were performed on first night in the K-prime filter {K'; Aq = 2.124/im; SX — 0.351/im), 
and on the second night in the J filter (Aq = 1.248/im; SX = 0.163;um). A single K' data 
frame was taken with an integration time of 2 seconds and 10 coadds; 10 frames were ac- 
quired for a total integration time of 200 seconds. A single J data frame was taken with an 
integration of 0.18s and 20 coadds; 10 frames were acquired for a total integration time of 
36 seconds. The weather on the night of the observations was poor with occasional heavy 
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clouds. 

The individual frames were background-subtracted and flat-fielded into a single final 
image for each filter. The central core of the resulting point spread functions had a width 
of FWHM = 0.07" 7 pixels) at J and FWHM = 0.09" 9 pixels) at K'. The final 
coadded images are shown in Figure [7} A faint source is detected 0.75" from the primary 
target at a position angle of PA = 129° east of north. The source is fainter than the primary 
target by A J = 4.70 ±0.05 mag and AK' = 3.95 ±0.05 mag. No other sources were detected 
within 5" of the primary target. 

If the faint companion is a dwarf star, the J — Ks color implies that it is a very late 
M-dwarf (~M5-M8; see Leggett et al. 2002; Ciardi et al. 2011) and would be at a distance 
of approximately 15±8 pc. If the companion is a giant star, the J — Ks color implies that it 
is an MO giant and would have an approximate distance of 10 kpc. Appendix A discusses in 
detail our use of the near-IR AO observation to convert the companion's brightness into Kp = 
14.4±0.2. The maximum line of sight extinction to the faint companion (as determined from 
the IRAS/DIRBE dust maps; see Schlegel et al., 1998) is Ay ~ 0.5 mag, which corresponds 
to an E{J — Ks) ~ 0.09 mag; such an excess would only change the implied spectral type 
by a single subclass (M5 dwarf and K5 giant) and would not appreciably change the derived 
distances. The red dwarf can be made much earlier in spectral type (and thus located at the 
same distance of HD 179070) but this would require significant reddening along the exact 
line of sight to the star which is ruled out based on the color excess listed for HD 179070 in 
the literature (E(b-y)=0.011). The primary target has a Hipparcos distance of 108 ± 10 pc; 
thus, the faint companion, whether it is a dwarf or giant, is not physically associated with 
the primary target. 

Source detection completeness was estimated by randomly inserted fake sources of vari- 
ous magnitudes in steps of 0.5 mag and at varying distances in steps of 1.0 FWHM from the 
primary target. Identification of sources was performed both automatically with the IDL 
version of DAOPhot and by eye. Magnitude detection limits were set when a source was not 
detected by the automated FIND routine or was not detected by eye. Within a distance of 
1 — 2 FWHM, the automated finding routine often failed even though the eye could discern 
two sources, particularly since the observations were taken in poor weather conditions. A 
summary of the detection efficiency as a function of distance from the primary star is given 
Table [Tj Beyond ~ 0.7", the detection limit is ~ 6 magnitudes fainter than the target. 
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4. Spectroscopic Observations 

Optical spectroscopy for HD 179070 was obtained by tliree ground-based telescopes as 
part of tlie Kepler mission follow-up program (Batalha et al. 2010). These observations 
include early reconnaissance spectra to assess stellar identification, rotation, and to make 
a first check on binary as well as detailed high precision radial velocity work using Keck 
HIRES. We compare all of our ground-based spectroscopic determined values in Table 2. 

Reconnaissance spectra of HD 179070 were obtained on 13 & 16 Sept. 2010 UT at the 
Kitt Peak 4-m telescope using the RCSpec instrument. The 4-m RCSpec setup used a 632 
1/mm grating (KPC-22b in second order) with a 1 arcsec slit to provide a mean spectral 
resolution of 1.6A per resolution element across the full wavelength range of 3750-5100A. 
The spectra were reduced in the normal manner with observations of calibration lamps and 
spectrophotometric stars (obtained before and after each sequence) and bias and flat frames 
collected each afternoon. Each fully reduced 4-m spectrum (see Figure 8) was cross-correlated 
and fit to both the entire MK standard stars digitally available in the "Jacoby Atlas" 
(Jacoby et al. 1984; covers all spectral and luminosity types) as well as to a suite of stellar 
models (ranging in Teg from 3500K to 7000K, log g from 1.0 to 5.0, and solar metallicity) 
available through the Spanish Virtual Observatorjj^ 

Spectral type, luminosity class, and other stellar parameters were provided by the best 
fit match and both of the 4-m spectra gave consistent results: F4-6 IV star with Teg = 
6250±250K, log g = 4.0±0.25, and metal poor ([Fe/H] = -0.15). No relevant vsini informa- 
tion was available from the moderate resolution (R~5000) 4-m spectra. 

As is common procedure for the Kepler Mission, all exoplanet candidate stars also 
receive high resolution, low signal to noise (S/N), spectroscopic observations to identify easily 
recognizable astrophysical false positives. One or two correctly timed spectra can help rule 
out many types of false positives, including single- and double-lined binaries, certain types of 
hierarchical triples and even some background eclipsing binaries, all of which show velocity 
variations and/or composite spectra that are readily detectable by the modest facilities used 
for these reconnaissance observations. We also use these spectra to estimate the effective 
temperature, surface gravity, metallicity, and rotational and radial velocities of the host star. 
Below is a brief description of the instrument, the data reduction, and the analysis performed 
in this step. 

We used the Tillinghast Reflector Echelle Spectrograph (TRES; Furesz 2008) on the 
1.5m Tillinghast Reflector at the Fred L. Whipple Observatory on Mt. Hopkins, AZ to 



■^httpt/Zsvo.laeff.inta.es/ 
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obtain a high resolution, low S/N spectrum of HD 179070 (S/N~7 per resolution element, 
360 second exposure) on 28 Sep 2010 UT. The observation was taken with the medium fiber 
on TRES, which has a resolving power of A/dA ~ 44,000 and a wavelength coverage of 3900- 
8900 angstroms. The spectrum was extracted and analyzed according to the procedures 
outlined by Buchhave et al. (2010). Cross-correlations were performed against the grid of 
CfA synthetic spectra, which are based on Kurucz models calculated by John Laird and rely 
on a hne-list compiled by Jon Morse. The template with the highest correlation coefficient 
yields an estimate of the stellar parameters: Tefr = 6250±125 K, log(g) — 4.0±0.25, and 
Vrot = 8± 1 km/s. The errors correspond to half of the grid spacing, although they neglect 
possible systematics, e.g., those introduced in the event the metallicity differs from solar. 
We find the absolute radial velocity to be Yrad — -19.1±0.3 km/s. 

Finally, HD 179070 was subjected to RV measurements with Keck (Vogt et al. 1994, 
Marcy et al. 2008). We note that the Ca-II K-line shows virtually no chromospheric reversal, 
giving an S value of 0.14, placing HD 179070 among the quietest G stars. We have obtained 
14 RV measurements with the Keck telescope and HIRES spectrometer (R=60,000). The 
exposure time for each spectrum was established by use of an exposure meter such that all 
would yield a consistent S/N and thus very similar RV precision. The Keck HIRES exposures 
were 150±30 sec in all cases. Each observation consisted of a triplet of exposures with each 
exposure having a signal-to-noise of 210 per pixel and an internal RV error of ~2 m s~^. We 
determined the RV from each exposure and took the weighted mean as the final RV for each 
triplet. The internal error (a typical uncertainty for each exposure) was based on the iodine 
line fits for each exposure. 

We adopted a jitter of 5 m s~^, typical for stars of such spectral type and rotational 
Vsini, adding the jitter in quadrature to the internal errors that were ~2 m s~^. Low-gravity 
F6 stars are well known to exhibit jitter of ~5-10 m/sec (Isaacson & Fischer 2010) due 
presumably to photospheric velocity fields but the exact origin remains unclear. F-type 
main sequence stars engage in a quasi-stable 5 Scuti-like phenomenon, a feature not seen 
in the G and K stars. Table 3 gives the times and the velocity measurements and their 
uncertainties. The jitter was not added to these uncertainties in quadrature, offering the 
reader a chance to see the uncertainties pre-jitter. The resulting 14 RVs had a standard 
deviation of 5.6 m s~^, consistent with the expected total errors. 

In Figure 9, we show the velocities of HD 179090 measured with the HIRES spectrometer 
on the Keck 1 telescope during ~90 days in 2010 and 2011. Sec Jenkins et al. (2011) and 
Batalha et al. (2011) for a detailed explanation of the (standard) method we used with the 
iodine cell to make these Doppler measurements. The velocities in Figure 9 present no clear 
long term variability on time scales of weeks or months. A periodogram of the velocities 
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reveals no significant power at any period from 0.5 d to tlie duration of tlie RV observations, 
85 dy, including periodicity in the velocities at the transit period of 2.7857 d. The Keck RVs 
are shown in Figure 10 as a function of known orbital phase (from the Kepler photometric 
light curve). The RV variation measured shows no modulation coherent with the orbital 
phase and is consistent with no change at all within the uncertainties. 

We performed a standard LTE spectroscopic analysis (Valenti and Piskunov 1996; 
Valenti and Fischer 2005) of a high resolution template spectrum from Keck-HIRES to 
derive an effective temperature, Tes = 6131 K, surface gravity, logg = 3.9 (cgs), metallicity, 
[Fe/H]=-0.05 and Vsini = 7.5 km/sec. Ground-based high spectral resolution support obser- 
vations of HD 179070 were also performed by the Kepler Asteroseismic Science Consortium 
(KASC). Two teams obtained estimates for [Fe/H]; Molenda-Zakowicz et al. (2010) listed 
two values (-0.15 and -0.23, and [Fe/H] = -0.15 was determined by H. Bruntt (priv. comm.) 
using NARVAL at Pic du Midi and the VWA analysis package. We adopt [Fe/H] =-0.15 for 
HD 179070 in this paper using the more reliable value recently obtained by Bruntt. 



5. Asteroseismic analysis 

5.1. Estimation of asteroseismic parameters 

HD 179070 was observed for one month by Kepler at a short cadence of 58.85 s. A 
time series was prepared for asteroseismic analysis in the manner described by Garcia et 



al. (2011). Fig. 11 plots the frequency-power spectrum of the prepared time series, which 
shows a beautiful pattern of peaks due to solar-like oscillations that are acoustic (pressure, 
or p) modes of high radial order, n. The observed power in the oscillations is modulated 
in frequency by an envelope that has an approximately Gaussian shape. The frequency of 
maximum oscillation power, z/max, has been shown to scale to good approximation as gT^^^'^ 
(Brown et al. 1991; Kjeldsen & Bedding 1995), where g is the surface gravity and T^s is the 
effective temperature of the star. The / identification, from visual inspection of the mode 
pattern in the frequency-power spectrum, is unambiguous and the n identification followed 
from the best-fitting to stellar evolutionary models (see below). The most obvious spacings in 
the spectrum are the large frequency separations, Az/, between consecutive overtones n of the 
same spherical angular degree, /. These large separations scale to very good approximation 
as (p)^''^, (p) oc M/R^ being the mean density of the star, with mass M and surface radius 
R (e.g. see Christensen-Dalsgaard 1993). 

Here, seven teams estimated the average large separation, (Az/), and z/max, using auto- 
mated analysis tools that have been developed, and extensively tested (e.g., see Campante 
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et al. 2010a; Christensen-Dalsgaard et al. 2010; Hekker et al. 2010; Huber et al. 2009; 
Karoff et al. 2010; Mosser & Appourchaux 2009; Mathur et al. 2010a; Roxburgh 2009) 
for application to the large ensemble of solar-like oscillators observed by Kepler (Chaplin et 
al. 2010,2011; Verner et al. 2011). A final value of each parameter was selected by taking 
the individual estimate that lay closest to the average over all teams. The uncertainty on 
the final value was given by adding (in quadrature) the uncertainty on the chosen estimate 
and the standard deviation over all teams. We add that there was excellent consistency 
between results, and no outlier rejection was required. The final values for (Au) and i^max 
were 60.86 ±0.55 yuHz and 1153±32/iHz, respectively. We did not use the average frequency 
separation between the / = and / = 2 modes (often called the small separation) in the 
subsequent modeling because HD 179070 turns out to be a subgiant (as indicated by, for 
example, the size of (Az/)), and in this phase of evolution the parameter provides little in 
the way of additional constraints given the modest precision achievable in it from one month 
of data (see Metcalfe et al. 2010; White et al., 2011). 

Use of individual frequencies increases the information content provided by the seismic 
data for making inference on the stellar properties. Six teams provided estimates of indi- 
vidual frequencies, applying "peak bagging" techniques developed for application to CoRoT 
(Appourchaux et al. 2008) and Kepler data (e.g., see Metcalfe et al. 2010; Campante et al. 
2010b; Mathur et al. 2010b; Fletcher et al. 2010). We implemented the procedure outlined 
in Campante et al. (2011) and Mathur et al. (2011) to select from the six sets of estimated 
frequencies one set that would be used to model the star. This so-called "minimal frequency 
set" contains estimates on modes for which a majority of the teams' estimates were retained 
after applying Peirce's criterion (Peirce 1852; Gould 1855) for outlier rejection. Use of one of 
the individual sets, as opposed to some average over all sets, meant that the modeling could 
rely on an easily reproducible set of input frequencies. The selected frequencies are listed in 
Table H 



5.2. Estimation of stellar properties 

We adopted two approaches to estimate the fundamental properties of HD 179070. In 
the first we used a grid-based approach, in which properties were determined by searching 
among a grid of stellar evolutionary models to get a best fit for the input parameters, which 
were (Az/), z/max, and Teg = 6131 ±44K, from spectroscopic observations made on the Keck 
telescope in support of the HD 179070 analysis, and [Fe/H]= —0.15 ± 0.06 from Molenda- 
Zakowicz et al. (2011). Descriptions of the grid-based pipelines used in the analysis may be 
found in Stello et al. (2009), Basu et al. (2010), Quirion et al. (2010) and Gai et al. (2011). 
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In the second approach, the individual frequencies Uni were analyzed by the Asteroseis- 
mic Modeling Portal (AMP), a web-based tool tied to TeraGrid computing resources that 
uses the Aarhus stellar evolution code ASTEC (Christensen-Dalsgaard 2008a) and adiabatic 
pulsation code ADIPLS (Christensen-Dalsgaard 2008b) in conjunction with a parallel genetic 
algorithm (Metcalfe & Charbonneau 2003) to optimize the match to observational data (see 
Metcalfe et al. 2009, Woitaszek et al. 2009 for more details). 

Each model evaluation involves the computation of a stellar evolution track from the 
zero-age main sequence (ZAMS) through a mass-dependent number of internal time steps, 
terminating prior to the beginning of the red giant stage. Exploiting the fact that (Az/) 
is a monotonically decreasing function of age (see Metcalfe et al. 2009, and references 
therein), we optimize the asteroseismic age along each evolution track using a binary decision 
tree. The frequencies of the resulting model are then corrected for surface effects following 
the prescription of Kjeldsen et al. (2008). A separate value of is calculated for the 
asteroseismic and spectroscopic constraints, and these values are averaged for the final quality 
metric to provide more equal weight to the two types of observables. The optimal model is 
then subjected to a local analysis employing a modified Levenberg-Marquardt algorithm 
that uses singular value decomposition (SVD) to quantify the values, uncertainties and 
correlations of the final model parameters (see Creevey et al. 2007). 



5.3. Results on stellar properties 

Both approaches to estimation of the stellar properties yielded consistent results on the 
mass and radius of the star. The final estimates are M = 1.34 ± 0.01(stat)±0.06(sys) Mq 
and R = 1.86 ± 0.02(stat)±0.04(sys) R©. The statistical uncertainties come from the SVD 
analysis of the best-fitting solution to the individual frequencies. The spreads in the grid- 
pipeline results - which reflect differences in, for example, the evolutionary models and input 
physics - were used to estimate the systematic uncertainties. 

The grid pipelines, which used only the average seismic parameters, showed (in some 
cases) two possible solutions for the age of the star (one around 3 Gyr and another around 
4 Gyr or higher) . Use of the individual frequencies resolved this ambiguity by giving a best- 
fitting solution that clearly favored the younger model, the best estimate of the age being 
r = 2.84 ± 0.10(stat)±0.33(sys) Gyr. Fig. 12 shows good agreement between the frequencies 



of the best-fitting model and the observed frequencies. This echelle diagram (e.g., see Grec et 
al. 1983) plots the frequencies against those frequencies modulo the average large frequency 
separation of 60.86 /xHz. Overtones of the same spherical degree / are seen to align in near- 
vertical ridges. The observed frequencies from Table |4] are plotted in black, along with their 
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associated la uncertainties; while the best-fitting model frequencies are plotted in color. 

6. Kepler Photometry Transit Fits 

The raw Q1-Q5 light curve is presented in the top panel of Figure 1. The trends observed 
on various timescales are a combinination of astrophysical phenomenon and instrumental ar- 
tifacts. These trends were removed from the lightcurve as our lightcurve model did not 
account for such effects. The photometric time series was prepared for modeling by inde- 
pendently detrending each quarterly Kepler time-series. A cubic polynomial was fitted and 
removed and then filtered with a 5-day running median. Any observations that occurred 
during transit where masked out during the calculation of the polynomial fits or medians. 
Our model fits for the physical and orbital parameters of the planetary system. The tran- 
sit shape was described by the analytic formulae of Mandel & Agol 2002. We adopted a 
non-linear limb-darkening law (Claret). Coefficients were calculated by convolving Atlas-9 
spectral models with the Kepler bandpass using the adopted estimate of Tes from Table 2 
and the asteroseismic value of log g. Limb darkening coefficients are held fixed for all transit 
fits. We assumed a Keplerian orbit for the planet with zero eccentricity. Our model fitted for 
the period (P), epoch (TO), impact parameter (b), the mean stellar density (p*), the ratio of 
the planet and star radii (Rp/i?^), radial velocity amplitude (K) and the radial velocity zero 
point (7). A set of best fit model parameters was constructed by fixing p+ to the best value 
from asteroseimology. Minimization of Chi-squared was found using a Levenberg-Marquart 
method allowing P, TO, b, Rp/i?^, K and 7 

To estimate the error on each fitted parameter a hybrid-MCMC approach was used 
similar to Ford (2005). The asteroseismic value of p^, and its statical error was adopted as 
a prior. A Gaussian Gibbs sampler was used to identify new jump values of test parameter 
chains. The width of the Gaussian sampler was initial determined by the error estimates from 
the best-fit model. After 500-chains where generated the chain success rate was examined 
and the Gaussian width was rescaled using Equation 8 Gregory (2011). This process was 
repeated after the generation of each 500 chains until the success rate for each parameter 
was between 22 and 28%. At which point the Gaussian width was held fixed. 

To handle the large correlation between the model parameters a hybrid MCMC algo- 
rithm was adopted based on Gregory (2011). The routine works by randomly using a Gibbs 
sampler or a buffer of previously computed chain points to generate proposals to jump to 
a new location in the parameter space. The addition of the buffer allows for a calculation 
of vectorized jumps that allow for efficient sampling of highly correlated parameter space. 
After the widths Gibbs sampler stabilized, 200,000 chains where generated. The process was 
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repeated 3 additional times to test for convergence via a Gelman Rubin test (Gelman & 
Rubin 1992). 

The 4 chain sets where combined and used to calculate the median, standard deviation 
and la bounds of the parameter distribution centered on the median value. Adopting the 
asteroseismic errors on the stellar mass and radius we computed the planetary radius (Rp), 
inclination angle (i) and semi-major axis (a) and also the scaled semi-major axis (a/i?^) and 
transit-duration {Tdur) from the model parameter distributions and report all values in Table 
5. 



7. Addressing Blend Scenarios 

The lack of a clear Doppler detection needed for dynamical confirmation of the nature 
of the transit signals in HD 179070 requires us to address the possibility that they are the 
result of contamination of the light of the target by an eclipsing binary falling within the 
photometric aperture ("blend"). The eclipsing binary may be either in the background or 
foreground, or at the same distance as the target in a physically associated configuration 
(hierarchical triple). Furthermore, the object producing the eclipses may be either a star or 
a planet. 

We explore the wide variety of possible false positive scenarios using the BLENDER tech- 
nique (Torres et al. 2004, 2011; Freesin et al. 2011), which generates synthetic light curves 
for a large number of blend configurations and compares them with the Kepler photometry 
in a sense. The parameters considered for these blends include the masses (or spec- 
tral types) of the two eclipsing objects (or the size of the one producing the eclipses, if a 
planet), the relative distance between the binary and the target, the impact parameter, and 
the eccentricity and orientation of the orbit of the binary, which can affect the duration 
of the events. Our simulations explore broad ranges in each of these parameters, with the 
eccentricities for planetary orbits limited to the maximum value recorded for known tran- 
siting systems with periods as short as that of HD 179070 (we adopted a conservative limit 
of e < 0.4; see htt p : / / exoplanet . eu/) , and eccentricities for eclipsing binaries limited to 
e < 0.1 (Raghavan et al. 2010). Scenarios that give significantly worse fits than a true transit 
model fit (at the 3-a level) are considered to be rejected. While this rejection reduces the 
space of parameters for viable blends considerably, it does not eliminate all possible blends. 
Constraints from follow-up observations described previously (such as high-resolution imag- 
ing and spectroscopy) as well as multi-band photometry available for the target allow us 
to rule out additional areas of parameter space (see Figure 13). We then estimate the a 
priori likelihood of the remaining blends in the manner described in the next section. To 
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obtain a Bayesian estimate of the probability that the transit events are due to a bona-fide 
planet, we must compare the a priori likelihood of such a planet and of a false positive 
(odds ratio). We consider the candidate to be statistically "validated" if the likelihood of a 
planet is several orders of magnitude greater than that of a blendj^ For full details on the 
BLENDER procedure, we refer the reader to the references cited above. Examples of other 
Kepler candidates validated in this way include Kepler-9d (Torres et al. 2011), Kepler-lOc 
(Fressin et al. 2011), Kepler-llg (Lissauer et al. 2011), Kepler-18b (Cochran et al 2011), 
and Kepler-19b (Ballard et al. 2011). 



7.1. Background blends 

We examined first the case of background eclipsing binaries composed of two stars. 
Our detailed simulations with BLENDER indicate that false positives of this kind are not 
able to match the observed shape of the transit well enough (either in depth, duration of 
ingress /egress, or total duration), or else they feature significant ellipsoidal variations out of 
transit that are not seen in the Kepler photometry. The best-fitting scenario of this kind 
gives a match to the observations that is worse than that of a true transiting planet model 
at the 6a level, which we consider unacceptable. We also find that blends involving evolved 
stars (giants or subgiants) orbited by a smaller star are easily ruled out, as well as those 
with a main-sequence star eclipsed by a white dwarf. In both cases the companion induces 
strong curvature out of eclipse due to the short orbital period (2.78 days), and for giants the 
large stellar radius additionally requires a grazing "V" -shaped transit to match the observed 
duration. 

When the object producing the eclipses is a planet rather than a star, ellipsoidal vari- 
ations are negligible, and the shape of the echpses (further attenuated by the light of the 
target) can more easily match the observed shape for a large range of properties of the stars 
and planets involved. An illustration of the constraints provided by BLENDER for false pos- 



itives of this kind is shown in Figure 14 Following the BLENDER nomenclature we refer to 
the target star as the "primary", and to the components of the eclipsing pair as the "sec- 
ondary" and "tertiary" (in this case a planet). The figure shows the landscape (goodness 
of fit compared to a true transiting planet model) projected onto two of the dimensions of 
parameter space, corresponding to the mass of the secondary on the horizontal axis and the 



^In the context of this paper we reserve the term "confirmation" for the unambiguous detection of the 
gravitational influence of the planet on its host star (e.g., the Doppler signal) to establish the planetary 
nature of the candidate; when this is not possible, as in the present case, we speak of "validation" , which 
involves an estimate of the false alarm probability. 
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relative distance between the primary and the binary on the vertical axis. The latter is cast 
for convenience here in terms of the difference in distance modulus. The colored regions 
represent contours of equal goodness of fit, with the 3-0" contour indicated in white. Blends 
inside this contour give acceptable fits to the Kepler photometry, and are considered viable. 
They involve stars that can be up to 7 magnitudes fainter than the target in the Kepler 
passband (as indicated by the dashed green line in the figure, corresponding to background 
stars with AKp = 7), and that are transited by a planet of the right size to produce the 
measured signal. Also indicated are other constraints that rule out portions of parameter 
space otherwise allowed by BLENDER. The blue hatched region represents blends that have 
overall colors for the combined light as predicted by BLENDER that are either too red (left) or 
too blue (right edge) compared to the measured color of the target {r — Kg = 1.314 ± 0.035, 
adopted from the KIC; Brown et al. 2011), at the 3-a level. The green hatched area rep- 
resents blends that are bright enough (up to 4 mag fainter than the target) to have been 
detected in our high-resolution spectroscopy as a second set of lines (see Sect. 4). With these 
observational constraints the pool of false positives of this kind is significantly reduced, but 
many remain. We describe in Sect. 7.3 how we assess their frequency. 



7.2. Blends involving physically associated stars 

Hierarchical triple configurations in which the eclipsing object (tertiary) is a star are 
easily ruled out by BLENDER, as these configurations invariably lead to the wrong shape for 
a transit. However, stars physically associated with the target that are orbited by a planet 
of the appropriate size can still mimic the light curve well when accounting for dilution from 



the brighter star HD 179070. The map for this type of blend is seen in Figure 15 In this 
case the color of the blend is not a strong discriminant, as all of these false positives are 
predicted to have r—Kg indices similar to that of the target itself. The expected brightness of 
the companion stars, though, is such that most would have been detected spectroscopically 
[AKp < 4; green hatched exclusion region), unless their RV compared to the target is small 
enough that their spectral lines are blended with those of the main star. Based on our 
spectroscopic observations we estimate conservatively that we would miss such companions 
if they had radial velocities within ~15kms~^ of the RV of the target. These blends are not 
eliminated by any other observational constraint; we estimate their frequency below. 
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7.3. Validation of Kepler-21 b 

With the constraints on false positives afforded by the combination of BLENDER and 
other follow-up observations, we may estimate the a priori likelihood of a blend following a 
procedure analogous to that explained by Fressin et al. (2001). 

For blends involving background stars transited by a planet, this frequency will depend 
on the density of background objects near the target, the area around the target within which 
such stars would go undetected, and the rate of occurrence of planets of the appropriate 
size transiting those stars. We perform these calculations in half-magnitude bins, with the 
following ingredients: a) the Galactic structure models of Robin et al. (2003) to estimate the 
number density of stars per square degree, subject to the mass limits allowed by BLENDER; b) 
results from our adaptive optics observations to estimate the maximum angular separation 
(Pmax) at which companions would be missed, as a function of magnitude difference relative 
to the target {Kp = 8.224) properly converted to the Kp band, as described in the Appendix; 
and c) the overall frequency of suitable transiting planets that can mimic the signal. The 
size range for these planets, as determined in our BLENDER simulations, is 0.38-2.0 -Rjup. To 
estimate the frequency of such planets we make use of the list of 1235 planet candidates 
released by the Kepler Mission (Borucki et al. (2011), based on the first four months of 
observation by the spacecraft. While these objects have not all yet been confirmed because 
follow-up is still in progress, the false positive rate is expected to be relatively low (typically 
less than 10%; see l25|l and will not affect our results significantly. We therefore assume that 
all of them represent true planets, and that the census of Borucki et al. (2011) is complete 
for objects of this size (see below). The estimated frequency of these planets in the allowed 
radius range is /pianet ~ 0.19%. 

The results of our calculation for the frequency of blends involving background stars is 
presented in Table [6j Columns 1 and 2 give the magnitude range for background stars and 
the magnitude difference compared to the target; column 3 lists the range of allowed masses 



for the stars, based on our BLENDER simulations (see Figure 14); columns 4 and 5 list the 
mean star densities and pmax, respectively, and column 6 gives the number of background 
stars we cannot detect, and is the result of multiplying column 4 by the area implied by 
Pmax- Finally, the product of column 6 and the transiting planet frequency of 0.19% leads 
to the blend frequencies in column 7. The sum of these frequencies is given at the bottom 
under "Totals", and is 8.0 x 10"^. 

For blends involving physically associated stars with RVs within 15kms~^ of the target, 
which would go unnoticed in our spectroscopic observations, we estimate the frequency 
through a Monte Carlo experiment. We simulate companion stars in randomly chosen orbits 
around the target, and randomly assign them transiting planets in the appropriate radius 
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range as a function of the secondary mass (see Figure 15), according to their estimated 
frequencies from Borucki et al. (2011). We then determine what fraction of these stars 
would be missed because of projected angular separations below the (/.'OS detection threshold 
from our speckle observations, velocity differences relative to the target under 15kms~^, or 
because they would induce a drift in the RV of the target that is undetectable in our Keck 
observations (i.e., smaller than ±10 ms^^ over a period of 82 days; see Sect. 4). Binary 
orbital periods, eccentricities, and mass ratios were drawn randomly from the distributions 
presented by Raghavan et al. (2010), and the mass ratios used in combination with our 
estimate of the mass of HD 179070 to infer the mass of the physical companions. We adopt 
an overall binary frequency of 34% from the same source. 

Based on these simulations we obtain a frequency for this type of false positive of 1.17 x 
10~^. However, as seen in Figure [is} planets involved in blends with physically associated 
stars can be considerably smaller (~0. 12-0. 18 -Rjup) than those involved in background blends 
(0.38-2.0 i?jup), so we must consider the potential incompleteness of the census of Borucki 
et al. (2011) at the smaller planet sizes. To estimate this we performed Monte Carlo 
simulations in which we calculated the signal-to-noise ratio for each of the Kepler targets 
that would be produced by a central transit of a planet with the period of HD 179070 and 
with a given radius in the allowed range. Adopting the Kepler detection threshold of 7 for 
the signal-to-noise ratio Jenkins et al. (2010), we determined the fraction of stars for which 
such a planet could have been detected during the four months in which that sample was 
observed. We have assumed that the signal-to-noise ratio increases with the square root 
of the transit duration and with the square root of the number of transits, and that the 
data were taken in a continuous fashion (except for gaps between quarters). In this way 
we obtained a completeness fraction of about 65%, although this may be slightly optimistic 
given that some transits could have been missed due to additional interruptions in the data 
flow for attitude corrections and safe mode events. This brings the frequency of hierarchical 
triple blends to 1.8 x 10"®. 

The total blend frequency is then the sum of the two contributions (background stars and 
physically associated stars with transiting planets), which is 8x 10~^ + 1.8 x 10~® = 2.6 x 10~®. 

Finally, following the Bayesian approach outlined earlier, we require also an estimate 
of the likelihood of a true planet around HD 179070 ("planet prior") to assess whether 
it is sufficiently larger than the likelihood of a blend, in order to validate the candidate. 
To estimate the planet prior we may appeal once again to the catalog of 1235 candidates 
from Borucki et al. (2011), which contains 99 systems with planetary radii within 3a of 
the measured value for HD 179070 {Rp = 1.64 ± 0.04 R^). The 3-a limit used here is for 
consistency with a similar criterion adopted above in BLENDER. Given the total number of 
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156,453 Kepler targets from which the 1235 candidates were drawn, we obtain a planet 
frequency of 99/156,453 = 6.3 x 10~^. Applying the same incompleteness factor described 
above, which holds also for the radius of this candidate, we arrive at a corrected planet prior 
of 9.7 X 10~^. This conservative figure is ~370 times larger than the blend frequency, which 
we consider sufficient to vahdate the planet around HD 179070 to a high degree of confidence. 
We note that this odds ratio is a lower limit, as we have been conservative in several of our 
assumptions. In particular, for computing the frequency of planets transiting background 
stars (Table [6]) we have included objects with sizes anywhere between the minimum and 
maximum planet radius allowed by BLENDER for stars of all spectral types (0.38-2.0 i?jup), 
whereas the planet size range for secondaries of a given mass is considerably smaller. This 
would reduce the frequency of this type of false positive, strengthening our conclusion. 



8. Limits to the Density and Mass for Kepler-21b 

To determine a statistically firm upper limit to the planet mass, we carried out an 
MCMC analysis of the Keck radial velocities with a Keplerian model for the planet's orbit. 
The resulting 2-sigma upper limit to the mass yielded the Keplerian model shown in Figures 
9 & 10 and gave the following upper limits: RV amplitude of K <3.9 m s~^; a planet mass 
of M = 10.4 Msarth (2cr), and a corresponding density of p < 12.9 g cm^'^. This upper 
limit to the density of 12.9gcm~^ is so high that the planet could be (compressed) solid or 
composed of admixtures of rocky, water, and gas in various amounts, unconstrained by this 
large upper limit to density. The 1-sigma upper limit to density is 7Agcm~^, still consistent 
with all types of interior compositions and would yield a planet mass of ~5.9 MEarth- 

If Kepler-21b contains a large rocky core, the high pressure inside such a massive planet 
would cause the silicate mantle minerals to compress to dense phases of post-perovskite; 
the iron core is also at higher density than inside Earth (Valencia et al. 2007). However, 
Kepler-21b could also have a small rocky core, be mostly gas and not be nearly as massive. 
The maximum core fraction expected for rocky planets of this radius corresponds to a planet 
with mass of 10.0 MEarth and a mean density of 12.5 g/cc (see mantle stripping simulations 
by Marcus et al. 2010) with a corresponding RV semi-amplitude of 2.3 m/sec, still below 
our detection limit. If Kepler-21b is a water planet with low silicate-to-iron ratio and 50% 
water by mass, its mass would be merely 2.2 MEarth, similar to that of Kepler-llf, but at 
mean density of 2.7 g/cc. The measured radial velocities provide neither a confirmation nor 
a robust limit (~10 Earth-masses) on the mass of Kepler-21b but suggests an upper limit 
near that of the maximum rocky core fraction theoretically allowed. The radial velocities 
certainly rule out higher mass companions (additional planets or stellar companions) with 
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orbital periods in the period range of up to approximately 200 days as any RV trend caused 
by such a companion would be apparent in Figs. 9 & 10. 

9. Conclusion 

Kepler photometry of the bright star HD 179070 reveals a small periodic transit-like 
signal consistent with a 1.6 REarth exoplanet. The transit signal repeats every 2.8 days and 
the complete phased light curve shows all of the events to be consistent in phase, amplitude, 
and duration. Analysis of the Kepler image data and difference images are well matched by 
model fits. Detailed point response function (PRF) models conclude that the source of the 
transit event is centered on or near to the center of HD 179070 itself. Furthermore, these 
models show that many faint background eclipsing binary scenarios, capable of blending 
light with that from HD 179070 to produce the transit-like event, can be eliminated. 

High resolution ground-based optical speckle imaging reveals no nearby companion star 
to within 5 magnitudes of HD 179070 itself. Near-IR AO observation, however, reveal a faint 
companion star 0.75" away and ~4 magnitudes fainter in K. Using a color transformation, 
this star is expected to be R~14.2, just below the detection limit of the speckle results. Spec- 
troscopic observations also confirm that no bright star is present near HD 179070 (within 
0.5"). Asteroseismology was performed for HD 179070 using the Kepler light curves. Adopt- 
ing the spectroscopically determined values for Tefr, log g, and [Fe/H], the mass, radius, and 
age of HD 179070 were well determined. 

Putting all of the above observations and models together, we conclude that the cause of 
the periodic transit event is indeed a small 1.6 REarth exoplanet orbiting the subgiant star HD 
179070. Transit models were fit to the highly precise Kepler hght curve data revealing that 
the exoplanet orbits every 2.78 days at an inclination of 82.5 degrees. The exoplanet has an 
equilibrium temperature near 1900K and is located 0.04 AU from its host star. Kepler-21b 
has been validated by detailed modeling of blend scenarios as a true exoplanet at greater 
than 99.7% confidence. We can only determine an upper mass limit for the exoplanet, ~10 
MEarth, resulting in a upper limit to the mean density of ~13 g cm^^. 

if ep/er continues to monitor HD 179070 and will eventually build up higher S/N phased 
transit light curves. These long term observations may allow other planets within this 
same system may be directly detected or detected via transit timing variations. Given the 
brightness of HD 179070, it is likely that continued radial velocity monitoring will take 
place with Keck or other current or planned radial velocity instruments. Given a consistent 
level of instrumental precision, the observed stellar jitter will slowly be averaged out and 
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velocity signals from this or other planets orbiting the host star may be detected. Finally, 
using a technique such as described in Schuler et al. (2011) high-resolution, high signal-to- 
noise echelle spectroscopy will provide detailed metal abundance values of the host star's 
atmosphere which may hold clues as to the formation, or not, of planetary bodies. 

We thank John Johnson for use of some of his Keck Time. This research has made 
use of the NASA/IPAC/NExScI Star and Exoplanet Database, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with the National 
Aeronautics and Space Administration. The authors would like to thank the Kepler Science 
Office and the Science Operations Center personal for their dedicated effort to the mission 
and for providing us access to the science office data products. The ground-based observa- 
tions reported on herein were obtained at Kitt Peak National Observatory, National Optical 
Astronomy Observatory, which is operated by the Association of Universities for Research 
in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 
Kepler was selected as the 10th mission of the Discovery Program. Funding for this mission 
is provided by NASA. 



10. Appexdix A: Transformation of Infrared Colors 

To understand the contribution of the faint infrared companion to the light curve in 
the Kepler bandpass, we need to convert the measured infrared color [J — K') to a Kepler 
magnitude (Kp). To do this, we have derived a color-color relationship {Kp—Ks vs. J — Ks) 
utilizing the Kepler targets from Ql public release and the photometry from the Kepler Input 
Catalog (KIC; Brown et al. 2011). Separating the KIC into dwarfs and giants as described by 
Ciardi et al. (2011), we have fitted the color-color relationship with a 5th-order polynomial 



for the dwarfs and a 3rd-order polynomial for the giants (see Figures 16 and 17) 



The dwarf and giant color-color relationships were determined separately from the 
Kepler magnitude and 2MASS magnitudes of 126092 dwarfs within the color range of 
—0.2 < J—Ks < 1.0 mag and 17129 giants within the color range of the —0.2 < J—Ks < 1.2 
mag. The resulting polynomial coefficients from the least squares fits for the dwarfs and gi- 
ants, respectively, are 



Dwarfs : Kp-Ks = 0.314377 + 3.85667a; + 3. 1761113;^ - 25.3126x^ + 40.7221x^ - 19.2112x^ 



Giants : Kp-Ks = 0.42443603 + 3.7937617a; - 2.3267277a;2 + 1.4602553a;^ 

where x = J — Ks. 
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The fits and the residuals are shown in Figures 16 and 17 the residuals for both the 
dwarfs and giants are well-characterized by gaussian distributions with means and widths of 
{Kpj_Ks - Kptrue) = -0.005 ± 0.083 mag and {Kpj_Ks - Kptme) = -0.002 ± 0.065 mag 
for the dwarfs and giants, respectively. The uncertainties in the derived Kepler magnitudes 
(Kp) are dominated by the physical widths of the color-color relationships. 

The real apparent photometry of the infrared companion was determined from the 
2MASS photometry of the primary target (J = 7.229 ±0.032 mag, Ks = 6.945 ±0.018 mag) 
which is a blend of the two sources. The above color-color relationships were determined 
using the Ks filter, but the observations were taken in the K' filter which has a slightly 
shorter central wavelength (2.148 /xm vs. 2.124 /xm). Typically, the Ks and K' filters yield 
magnitudes which are within 0.02 — 0.03 magnitudes of each other and have zero-point flux 
densities within 2% (AB magnitudes = 1.86 and 1.84 for Ks and K' , respectively; Tokunaga 
& Vacca 2005). Given the quality of the weather and resulting photometry, the width of 
the color-color relationships, and lack of an if-band observation to aid in the transformation 
of the K' observations, we have equated K' to Ks in these calculations and propagated an 
additional uncertainty of 0.03 mag in the derivation of the J — Ks color and the Kepler 
magnitude (Kp). 

Deblending the infrared photometry results in the following infrared magnitudes for the 
primary target and faint companion of J = 7.24±0.07 and 11.94±0.07 mag,respectively, and 
Ks = 6.97 ± 0.07 and Ks = 10.92 ± 0.07 mag, where we have propagated the uncertainties 
of the faint companion on to the photometry of both stars. The infrared colors of the 
faint companion is J — Ks = 1.0 ± 0.1 mag, which corresponds to a Kp — Ks color of 
Kp — Ks = 3.54 ± 0.14 mag if the faint star is a dwarf and a color of Kp — Ks = 3.35 ± 0.14 
mag if the faint star is a giant. 

Applying the deblended Ks magnitude of the companion {Ks = 10.92 ± 0.07 mag), 
we derive a magnitude for the faint star in the Kepler bandpass for the dwarf- and giant- 
star relationships of Kp = 14.5 ± 0.2 mag and Kp = 14.3 ± 0.2 mag, respectively. The 
companion is fainter than the primary target, in the Kepler bandpass, by AKp = 6.3 mag 
if the star is dwarf and AKp = 6.1 mag if the star is a giant. Note the primary star 
dominates the photometry in the Kepler aperture; after deblending the Kepler magnitude 
of the star changes from Kp = 8.224 mag to Kp = 8.227 mag if the companion is a dwarf 
or to Kp = 8.228 mag if the companion is a giant. 

The above relationships only work if both J and Ks magnitudes are known, but often 
only one of the filters is available. Being able to convert a single J and Ks magnitude into 
an expected Kp is extremely useful - particularly, for determining sensitivity limits for the 
AO imaging. Towards this end we have utilized the Kepler Input Catalog to determine the 
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expected Kp — J and Kp — Ks colors for a given J or Ks magnitude. 



Histograms of the Kp — J and Kp — Ks colors are shown in Figure 18 where the median 
and mode of the color are marked; the spread in color, as described by the dispersion of the 
colors, is fairly large. The medians, modes (both of which are delineated in the histograms), 
and dispersions of the colors are Kp — J = 1.477, 1.275, and 0.626 mag and Kp — Ks = 2.139, 
1.775, and 0.803 mag. 

It is not unexpected that the measured median colors would be dependent upon the real 
apparent infrared magnitude; as the photometry becomes more sensitive to fainter and fainter 
sources, more intrinsically fainter (and redder) sources should contribute more significantly 
to the color distribution. To explore this effect, we have computed the median color {Kp — J 



and Kp — Ks) as a function of the real apparent infrared magnitude (see Fig 18). The 
dispersion per bin is fairly large (0.4 — 0.9 mag), but the colors show smooth systematic 
trends as a function of magnitude with a range of 0.5 — 0.6 mag. 

We have characterized these curves with 5th-order polynomials and have done a linear 
extrapolation for magnitudes fainter than the data range: 

Kp-J = -398.04666+149.08127J-21.952130j2+1.5968619J^-0.057478947J^+0.00082033223J^ 
for (10 < J < 16.7 mag) 

iTp - J = 0.1918 + 0.08156J 

for (J > 16.7 mag) and 

Kp-Ks = -643.05169+246.00603Ks-37.136501i^s^+2.7802622Ks^-0.10349091K/+0.0015364343Js:s^ 
for (10 < Ks < 15.4 mag) 

Kp- Ks = -2.7284 + 0.3311i^s 

for {Ks > 15.4 mag). 

The trends seen in the color vs. magnitude relationships are not unexpected. At the 
brighter magnitudes, the distribution of stars is dominated by infrared bright stars (i.e. 
giants) and thus, are dominated by relatively red stars. As the magnitude limit is increased, 
the dwarf stars begin to contribute to the color distribute starting with the bluer (more 
luminous) stars and the median colors become bluer. As the magnitude limits are pushed 
even further, the intrinsically fainter (i.e., red) dwarf stars begin to dominate the sample, 
and the median colors become increasingly red as the magnitude limit is increased. 
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Table 1. Approximate Radial Source Sensitivity 



Distance 


Distance 


A 1 
ZA J 


Distance 


lAPi. S 


VV rliVi J 


\") 


(mag) 


( ii\ 

\") 


(mag) 


1 


0.07 


1.0 


0.09 


1.5 


2 


0.14 


1.5 


0.18 


2.0 


3 


0.21 


2.0 


0.27 


2.5 


4 


0.28 


2.5 


0.36 


3.0 


5 


0.35 


3.0 


0.45 


3.5 


6 


0.42 


3.5 


0.54 


4.0 


7 


0.49 


4.0 


0.63 


4.5 


8 


0.56 


4.5 


0.72 


5.0 


9 


0.63 


5.0 


0.81 


5.5 


10 


0.70 


5.5 


0.90 


6.0 


11 


0.77 


6.0 


0.99 


6.5 



Table 2. Kepler image centroid offsets in pixels for HD 17907 



Quarter^ 


Row offset 


Column offset 


Offset distance 


1 


2.09 ± 0.02 


1.23 ±0.03 


2.43 ± 0.02 


3 


-0.44 ±0.03 


1.18 ±0.07 


1.26 ±0.07 


4 


0.18 ±0.07 


-0.35 ±0.07 


0.39 ±0.07 


5 


-0.30 ±0.02 


0.35 ± 0.03 


0.46 ± 0.03 



'^Quarter 1 was short (1 month) and provides less reliable 
ccntriod measurements. Quarter 2 data was not included 
due to its excess noise. 
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Table 2. Spectral Analysis of HD 179070 



Source 


Te// (K) 


log g (cgs) 


[Fe/H] 


Vsini (km/sec) 


N(2004)^ 


6137 




-0.015 




KPNO 4-m 


6250±250 


4.0±0.25 


-0.15±0.15 




TRES 


6250±125 


4.0±0.25 


0.0±0.25 


8.0±1.0 


Keck HIRES 


6131±44 


3.9±0.1 


-0.05±0.1 


7.5±1.0 


MZ(2010)'' 


6063±126 


4.04±0.07 


-0.23±0.09 


<5 


^^(2010)"^ 


6145±65 


4.15±0.10 


-0.15±0.06 


<5 


HB^ 






-0.15 




Adopted 


6131±44 


4.0±0.1 


-0.15±0.06 


7.75±1.0 



^Nordstrom et al., 2004 

''Molenda-Zakowicz et al., 2010, two solutions listed 
'^H. Bruntt, private communication 
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Table 3. Radial Velocities for HD 179070 



BJD 


RV 


Unc. 


-2440000 


(m/sec j 


(m/sec) 


±o^ou .uoo 


9 fil 


9 1 


1 "1439 941 


8 44 


9 9 




Q QQ 


9 9 


1 ^AAf] 771 


n 

u.oo 


9 9 




n 7'? 


9 n 


1 ^zLzLD 77^ 


(\ 1 8 
u . ±o 


9 9 


1 ^zLzLD QR7 


qq 


9 

Zi .0 




^ q^ 


9 

Zi .0 


1 ^AACi QQ1 

X (J ^ U . t? t? X 


n 84 


9 9 

Zi .Zi 


X(J'4:(J(J.OZ(U 


o.ux 


9 1 

Zi . X 


1 ^4^^ 

X(J'4:(JU.OZiO 


-7 "^n 


9 1 

£1 . X 


XtJ^(J(J.OOU 


9 74 


9 n 


1 ^4(S4 7RR 


14 q6 


1 8 

X .0 


1 ^4(S4 7Qn 


1 D 4^ 
xu 


1 R 
X .0 


1 ^4(S4 7Q'^ 


1 n 79 

XU. 1 z, 


1 R 
X .0 


XtjTiUtJ.OUU 


14 74 


1 8 

± .0 


lo^yjo .oyjcf 


4 19 

XZi 


1 q 




7 47 


1 9 




-D "16 


1 9 


1 "izlfifi 797 


-■^ Qfi 


9 n 


i54DD.7z9 


1 OT 

1.87 


z.O 


15467.848 


-9.35 


1.9 


15467.849 


-0.91 


2.0 


15467.851 


3.14 


2.1 


15468.714 


-1.64 


2.0 


15468.716 


2.15 


1.9 


15468.718 


1.55 


2.0 


15469.753 


-4.89 


1.7 


15469.755 


-8.50 


1.8 


15469.758 


1.13 


1.8 
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Table 3 — Continued 



BJD 


RV 


Unc. 


-2440000 


(m/sec) 


(m/sec) 


15471.847 


1.18 


2.0 


15471.849 


-3.06 


1.8 


15471.852 


-7.47 


1.8 


15486.819 


-2.46 


2.0 


15486.822 


0.56 


2.2 


15486.824 


-2.85 


2.0 


15490.819 


-3.67 


2.0 


15490.821 


-5.84 


2.1 


15490.823 


-7.85 


1.9 


15521.759 


-8.33 


2.2 
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Table 4. Estimated frequencies Uni of HD 179070 (in /xHz). 



n 


/ = 


1^1 


/ = 2 


12 






850.13 ±3.57 


13 


855.05 ± 2.09 


885.51 ± 1.23 


907.92 ± 4.93 


14 


918.26 ±1.16 


946.67 ± 1.28 


975.03 ± 1.31 


15 


979.45 ±0.25 


1005.25 ±0.75 


1034.88 ± 1.70 


16 


1039.33 ± 0.42 


1064.76 ±0.72 


1095.85 ± 1.23 


17 


1098.37 ±0.86 


1125.63 ±0.44 


1155.74 ± 1.46 


18 


1159.28 ±0.88 


1187.41 ±0.58 


1215.74 ± 1.95 


19 


1221.45 ± 1.03 


1248.53 ±0.61 


1279.21 ± 1.69 


20 


1282.74 ± 0.85 


1308.73 ± 0.71 


1339.38 ± 1.45 


21 


1341.48 ± 0.56 


1370.70 ± 1.02 


1399.88 ±2.19 


22 


1404.24 ± 1.85 


1432.05 ± 1.72 
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Table 5. Transit Model Parameters (e=0) 







A/To/H ion 
iVit^Cllclll 


oiciev 


4-1 rr 


1 n 


Adopted Values 












(Mo) 


1.340 




0.010 


0.010 


-0.010 




1.860 




0.020 


0.020 


-0.020 




4.0190 


4.0196 


0.0090 


0.0087 


-1.0694 


(g/cm^) 


0.2886 


0.2891 


0.0087 


0.0077 


-0.0102 


Derived Values 












Rp (Rj) 


0.1459 


0.1456 


0.0035 


0.0034 


-0.0038 


P (days) 


2.785755 


2.785755 


0.000032 


0.000031 


-0.000034 


i (deg) 


82.58 


82.59 


0.29 


0.28 


-0.31 


TO* 


193.8369 


193.8368 


0.0016 


0.0016 


-0.0016 




4.910 


4.913 


0.050 


0.043 


-0.058 




0.00806 


0.00804 


0.00018 


0.00018 


-0.00019 


b 


0.640 


0.639 


0.023 


0.020 


-0.028 


a (AU) 


0.042507 


0.042509 


0.000106 


0.000098 


-0.000119 


Trfur (h) 


3.438666 


3.438982 


0.078588 


0.070336 


-0.091437 


T 


1956±297 











and R^ are fixed to asteroseismic values 
tT0=BJD-2454900 
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Table 6. Blend frequency estimate for HD 179070 for scenarios involving background stars 

transited by a planet. 



Kp range 


AKp 


Stellar mass 


Stellar density 


pmax 


Stars 


Blends'^ 


(mag) 


(mag) 


range (A/-:-) 


(pci sc}. (leg) 


(") 


(xlQ-*') 


(xlQ-*') 


(1) 


(-') 


(•■-!) 


(1) 




((') 


(') 


8.2-8.7 


0.5 












8.7-9.2 


1.0 












9.2-9.7 


1.5 












9.7-10.2 


2.0 












10.2-10.7 


2.5 












10.7-11.2 


3.0 












11.2-11.7 


3.5 












11.7-12.2 


4.0 












12.2-12.7 


4.5 


0.80-1.40 


115 


0.50 


6.97 


0.0132 


12.7-13.2 


5.0 


0.80-1.40 


204 


0.60 


17.8 


0.0338 


13.2-13.7 


5.5 


0.80-1.40 


267 


0.75 


36.4 


0.0692 


13.7-14.2 


6.0 


0.80-1.40 


391 


0.85 


68.5 


0.130 


14.2-14.7 


6.5 


0.80-1.28 


438 


0.95 


95.8 


0.182 


14.7-15.2 


7.0 


0.85-1.15 


733 


1.05 


195.9 


0.372 


15.2-15.7 


7.5 












15.7-16.2 


8.0 












Totals 






2148 




421.4 


0.800 



Total blend frequency = 8.0 x 10~ 



*The range of radii allowed by BLENDER for the planets involved in these blends is 0.38— 
2.0i?jup, and the planet frequency used for the calculation is /planet = 0.19% (see text). 

Note. — Magnitude bins with no entries correspond to brightness ranges in which BLENDER 
excludes all blends, or that are ruled out by spectroscopic constraints. 
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Fig. 1. — Kepler light curve of HD 179070 covering Quarters to 5. The top raw hght 
curve covers 164 separate transit events for the small exoplanet orbiting the star. The 
middle panel shows a typical normalized section from the full light curve in which transits 
are visible (positions marked with dotted lines). The bottom panel shows the detrended, 
binned, and phase folded-data (see §6) overplotted by our model fit (red line). 
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Fig. 2. — A typical Kepler Quarter 5 pixel image of HD 179070 showing the saturation 
spilled along columns and the non-saturated wings around the core. 



difference image 



690 



720 




608 610 612 614 



Fig. 3. — Kepler Quarter 5 average difference image of HD 179070 showing the pixels that 
change during the transit event. 
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Fig. 4. — Left: Kepler Quarter 5 average difference image with the saturated pixels set to 
zero. Right: The corresponding PRF model difference image. 
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Fig. 5. — Quarter 5 model difference images in which the modeled faint background eclipsing 
binary star was offset in the column direction from the center location of HD 179070. As the 
column offset increases (left to right) the difference signal in the left most saturated pixels 
goes away. A similar result is seen for a column shift to the left. 
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Fig. 6. — Speckle observation of HD 179070 showing that no hne-of-sight or real companions 
exist from 0.05 to 2.8 arcsec of the star to a limit of 5.3 magnitudes in R (5.0 magnitudes 
in V) fainter than the star itself. The reconstructed images at the bottom of the plot have 
N up and E left. The horizontal line in the top plot shows the 5 sigma detection limit 
for companions against the sky background (open squares) and the vertical line at 0.2 arc 
seconds is added to show the inner limit for conservative multi-fringe speckle detections. 
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Fig. 7.— (top) ARIES AO images of HD 179070 in J (left) and Ks (right). The inner 2" 

_1 _1 JJ.1_ 1 J._ 1\T T71 rnl_ _ 1 _J J_ _1 1__ _ J._ J.1_ _ 



-44- 



1 ' 1 ' 1 1 1 ' r 




I I I I I \ I I I 1 I I 

3600 3800 4000 4200 4400 4600 4800 

Wavelength (A) 



Fig. 8. — The Kitt Peak 4-meter spectrum of HD 179070 obtained on 16 September 2010 
UT. The F6 IV star has an effective temperature of 6131K and a log g of 4.0. See Table 2. 
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Fig. 9. — Radial velocity measurements for HD179070 from the Keck HIRES spectrometer 
arc presented as a function of time. Internal errors of ~2m/s are added in quadrature 
to 5m/s of jitter to account for uncertainty in the measurements due to spectral type. A 
Keplerian orbit for Keplcr-21b is ovcrplotted (see test). The radial velocity amplitude does 
not correlate to the expected phase of the planet nor is any additional coherent variation 
observed. The radial velocities therefore provide only an upper limit to the amplitude. The 
small RMS scatter of 5 m s~^ imposes an upper limit on the planet mass (see text) and rules 
out a grazing incidence eclipsing binary. 
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Fig. 10. — Keck velocities vs orbital phase for HD 179070. The velocities exhibit no evidence 
of coherence with orbital phase, thus offering only a 3a upper limit to the mass of the planet 
of ~20 Earth-masses. The solid hne shows the expected RV curve for a 10.4 Earth-mass 
planet orbiting HD 179070 (see §8). 
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Fig. 11. — Frequency-power spectrum of HD 179070, showing a rich pattern of overtones of 
solar-hke oscillations. The rising background toward lower frequencies is due to convective 
granulation. 
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Fig. 12. — Echelle diagram showing the observed frequencies from Table |4| along with their 
la uncertainties (black symbols with error bars), and the best-fitting model frequencies 
(color). Different symbol styles denote different spherical degrees, /, circles showing / = 0, 
triangles / = 1 and squares 1 = 2. 
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Fig. 13. — A schematic that demonstrates various techniques that rule out the brightness of 
potential blends in the Kepler aperture. The yellow star in the upper left corner represents 
HD 179070. The red star represents Kepmag=14.5 the companion discovered at a separa- 
tion of 0.7" with AO. This star was just undetected with speckle imaging due to the late 
stellar type. Modeling of the transiting object rules out a blend by any star fainter than 
Kepmag=15.8 and a conservative estimate of the Kepler centroids eliminates any blend with 
a separation greater than 4" . 
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Fig. 14. — Map of the surface (goodness of fit) for blends involving background stars 
transited by a larger planet. The vertical axis represents the distance between the background 
star and the target HD 179070, expressed in terms of the difference in the distance modulus. 
Only blends inside the solid white contour match the Kepler light curve within acceptable 
limits (3(7, where sigma is the significance level of the difference compared to a transit 
model fit; see Fressin et al. 2011). Lighter-colored areas (red, orange, yellow) mark regions 
of parameter space giving increasingly worse fits to the data (4cr, 5cr, etc.), and correspond 
to blends we consider to be ruled out. The hatched blue region on the lower left corresponds 
to blends that can be excluded because of their overall r — Ks colors, which are too red 
compared to the measured index for HD 179070, by more than 3a (0.10 mag). A smaller 
similar region is visible on the right. Blends that are bright enough to have been detected 
spectroscopically are indicated by the hatched green area, corresponding to contaminating 
stars that are up to 3 mag fainter than the target. The faintest blends that remain can be 
as much as 7 mag dimmer than the target (dashed green line). 
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Fig. 15. — Similar to Figure 14 for the case of hierarchical triple systems in which the 
secondary star is transited by a planet, mimicking the signal in HD 179070. In this case 
the vertical axis represents the radius of those planets. Blends inside the white 3-cr contour 
have light curves that match the shape observed for HD 179070. While the r — Ks colors 
of these blends are indistinguishable from that of a single star like HD 179070 at the 3-cr 
level, the stars involved are all bright enough that most would have been detected in our 
high-resolution spectra as a second set of lines. This is indicated by the green hatched area. 
Only those with a RV such that the lines are completely blended with those of the target 
would escape notice (see text). 
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Fig. 16. — Top: Kepmag — Ks vs J — Ks for Ql dwarf stars. The grey line represents the 
5th order polynomial fit. Bottom: histogram of the residuals from the polynomial fit. 
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Fig. 17. — Top: Kepmag — Ks vs J — Ks for Ql giant stars. The grey line represents the 
3rd order polynomial fit. Bottom: histogram of the residuals from the polynomial fit. 



- 54 - 




Fig. 18. — Top: (left) Histogram of the Kp — J color for all sources in the KIC. The dashed 
line marks the median value and the dashed-dot line marks the peak of the histogram, 
(right) Median Kp — J color for all sources in the KIC as a function of J mag (binsize = 
0.02 mag in J). The horizontal dashed and dashed-dot lines mark the same values as in 
the histogram. Three representative dispersions per bin are shown as error bars. The white 
solid line represents a 5th-order polynomial fit to the data, and the dotted line represents a 
linear fit to sources fainter than J — Ih mag extrapolated to magnitudes beyond the limits 
of the KIC. Bottom: (left) Histogram of the Kp — Ks color for all sources in the KIC. (right) 
Median Kp — Ks color for all sources in the KIC as a function of Ks mag (binsize = 0.02 
mag in Ks). The various lines demark the same values as in the top plots but for Kp — Ks 
values. 



